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Matched-filtering for the identification of compact object mergers in gravitational-wave antenna 
data involves the comparison of the data stream to a bank of template gravitational waveforms. 
Typically the template bank is constructed from phenomenological waveform models since these 
can be evaluated for an arbitrary choice of physical parameters. Recently it has been proposed that 
singular value decomposition (SVD) can be used to reduce the number of templates required for de- 
tection. As we show here, another benefit of SVD is its removal of biases from the phenomenological 
templates along with a corresponding improvement in their ability to represent waveform signals 
obtained from numerical relativity (NR) simulations. Using these ideas, we present a method that 
calibrates a reduced SVD basis of phenomenological waveforms against NR waveforms in order to 
construct a new waveform approximant with improved accuracy and faithfulness compared to the 
original phenomenological model. The new waveform family is given numerically through the in- 
terpolation of the projection coefficients of NR waveforms expanded onto the reduced basis and 
provides a generalized scheme for enhancing phenomenological models. 

PACS numbers: 04.30.-w, 04.25.D-, 04.25.dg 



I. INTRODUCTION 

Developments are currently underway to promote the 
sensitivity of LIGO and to improve its prospect for de- 
tecting gravitational waves emitted by compact object 
binaries [D 12] • Of particular interest are the detection 
of gravitational waves released during the inspiral and 
merger of binary black hole (BBH) systems. Detection 
rates for BBH events are expected to be within 0.4-1000 
per year with Advanced LIGO [3]. It is important that 
rigorous detection algorithms be in place in order to max- 
imize the number of detections of gravitational wave sig- 
nals. 

The detection pipeline currently employed by LIGO 
involves a matched-filtering process whereby signals are 
compared to a pre-constructed template bank of gravi- 
tational waveforms. The templates are chosen to cover 
some interesting region of mass-spin parameter space and 
are placed throughout it in such a way that guarantees 
some minimal match between any arbitrary point in pa- 
rameter space and its closest neighbouring template. Un- 
fortunately, the template placement strategy generally 
requires many thousands of templates (e.g. [1]) evalu- 
ated at arbitrary mass and spin; something that cannot 
be achieved using the current set of numerical relativity 
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(NR) waveforms. 

To circumvent this issue, LIGO exploits the use of an- 
alytical waveform families like phenomenological models 
[6] or effective-one-body models [JJ [8]. We shall fo- 
cus here on the Phenomenological B (PhenomB) wave- 
forms developed by JB]. This waveform family describes 
BBH systems with varying masses and aligned-spin mag- 
nitudes (i.e. non-precessing binaries). The family was 
constructed by fitting a parameterized model to existing 
NR waveforms in order to generate a full inspiral-merger- 
ringdown (IMR) description as a function of mass and 
spin. The obvious appeal of the PhenomB family is that 
it allows for the inexpensive construction of gravitational 
waveforms at arbitrary points in parameter space and can 
thus be used to create arbitrarily dense template banks. 

To optimize computational efficiency of the detection 
process it is desirable to reduce the number of tem- 
plates under consideration. A variety of reduced bases 
techniques have been developed, either through singular- 
value decomposition (SVD) jHHU], or via a greedy algo- 
rithm [11] . SVD is an algebraic manipulation that trans- 
forms template waveforms into an orthonormal basis with 
a prescription that simultaneously filters out any redun- 
dancies existing within the original bank. As a result, 
the number of templates required for matched-filtering 
can be significantly reduced. In addition, it has been 
shown in |12j that, upon projecting template waveforms 
onto the orthonormal basis produced by the SVD, in- 
terpolating the projection coefficients provides accurate 
approximations of other IMR waveforms not included in 
the original template bank. 

In this paper, we continue to explore the use of the 
interpolation of projection coefficients. We take a novel 
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approach that utilizes both the analytic PhenomB wave- 
form family [6] and NR hybrid waveforms 13 15] . We ap- 
ply SVD to a template bank constructed from an analyt- 
ical waveform family to construct an orthonormal basis 
spanning the waveforms, then project the NR waveforms 
onto this basis and interpolate the projection coefficients 
to allow arbitrary waveforms to be constructed, thereby 
obtaining a new waveform approximant. We show here 
that this approach improves upon the accuracy of the 
original analytical waveform family. The original wave- 
form family shows mismatches with the NR waveforms as 
high as 0.1 when no extremization over physical parame- 
ters is applied (i.e., a measure of the "faithfulness" of the 
waveform approximant), and mismatches of 0.02 when 
maximized over total mass (i.e., a measure of the "effec- 
tualness" of the waveform approximant). With our SVD 
accuracy booster, we are able to construct a new wave- 
form family (given numerically) with mismatches < 0.005 
even without extremization over physical parameters. 

This paper is organized as follows. We begin in Sec- 
tion [TT] where we provide definitions to important ter- 
minology used in our paper. We then compare our NR 
hybrid waveforms to the PhenomB family and show that 
a mass-bias exists between the two. In Section IIIII we 
present our SVD accuracy booster applied to the case 
study of equal-mass, zero-spin binaries. In Section |TV| we 
investigate the feasibility of extending this approach to 
include unequal-mass binaries. We finish with concluding 
remarks in Section fVl 



II. GRAVITATIONAL WAVEFORMS 



complicate our approach by making waveforms longer in 
the low frequency domain. 

As a measure of the level of agreement between two 
waveforms, h(t) and g(i), we will use their match, or 
overlap, C(h,g) flBHTg] . We define 



0(h,g) 



max 

AT 



(h,g) 



|g| 



(2) 



where (h, g) is the standard complex inner product and 
the norm ||h|| = \J (h, h) . We always consider the over- 
lap maximized over time- and phase-shifts between the 
two waveforms. The time-maximization is indicated in 
([2]), and the phase-maximization is an automatic conse- 
quence of the modulus. Note that < C(h, g) < 1. For 
discrete sampling at points = to + iAt we have that 



*(ti) 



(3) 



where g*(t) is the complex conjugate of g(t). Without 
whitening, ^ would need to be evaluated in the fre- 
quency domain with a weighting factor 1/S n (f). The 
primary advantage of ^ is its compatibility with formal 
results for the SVD, which will allow us to make more 
precise statements below. When maximizing over time- 
shifts AT, we ordinarily consider discrete time-shifts in 
integer multiples of At, as this avoids interpolation. Af- 
ter the overlap has been maximized, it is useful to speak 
in terms of the mismatch, A^(h, g), defined simply as 



A4(h,g) = l-0(h,g). 



(4) 



A. Terminology 



We use this quantity throughout the paper to measure 
the level of disagreement between waveforms. 



A gravitational waveform is described through a com- 
plex function, h(f), where real and imaginary parts store 
the sine and cosine components of the wave. The specific 
form of h(t) depends on the parameters of the system, 
in our case the total mass M = mi + 7712 and the mass- 
ratio q = 7711/7712. While h(i) is a continuous function of 
time, we discretize by sampling h(i^), where the sampling 
times ti have uniform spacing At — 2~ 15 s. 

We shall also whiten any gravitational waveform h(t). 
This processes is carried out in frequency space via 

where S n (f) is the LIGO noise curve and h(/) is the 
Fourier transform of h(t). The whitened time-domain 
waveform, h w (f), is obtained by taking the inverse 
Fourier transform of ([!]). In the remainder of the paper, 
we shall always refer to whitened waveforms, dropping 
the subscript "w". For our purposes it suffices to take 
S n (f) to be the Initial LIGO noise curve. Using the Ad- 
vanced LIGO noise curve would only serve to needlessly 



B. NR Hybrid Waveforms 

We use numerical waveforms computed with the Spec- 
tral Einstein Code SpEC [19] . Primarily, we use the 15- 
orbit equal- mass (mass-ratio q = 1), zero-spin (effective 
spin x = 0) waveform described in [131 120] . In Section 
|IV[ we also use unequal mass waveforms computed by 
|14) . The waveforms are hybridized with a TaylorT3 
post-Newtonian (PN) waveform as described in [T5J [ST] 
at matching frequencies Mui = 0.038,0.038,0.042,0.044, 
and 0.042 for mass-ratios q = 1,2,3,4, and 6, respec- 
tively. 

TaylorT4 at 3.5PN order is known to match NR sim- 
ulations exceedingly well for equal-mass, zero-spin BBH 
systems [20] (see also Fig. 9 of 15J ) . For q= 1, a Tay- 
lorT3 hybrid is very similar to a Taylor T4 hybrid, cf. 
Figure 12 of [T3]. The mismatch between TaylorT3 and 
TaylorT4 hybrids is below 10~ 3 at M = 10 M , drop- 
ping to below 10~ 4 for 15 M < M < 20 M , and 10~ 5 
for 20 M < M < 100 M . These mismatches are sig- 
nificantly smaller than mismatches arising in the study 
presented here, so we conclude that our results are not 
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FIG. 1: The dashed line in the top panel traces the mis- 
match between equal-mass, zero-spin NR and PhenomB wave- 
forms of the same total mass M. Mismatch is reduced 
(solid line) by searching to find the mass M' for which 
.M[hNR.(M), Iipb(M')] is a minimum. We generally find that 
M' < M as shown in the bottom panel where the solid line 
traces the mass-bias (M — M')/M. For comparison, the dot- 
dashed curve in the bottom panel traces the mass spacing 
(Mk — Mfe_i)/Mfc for a template bank of PhenomB waveforms 
satisfying 7W[h PB (M fc _i), h PB (M fe )] = 8~ 2 /2. 



influenced by the accuracy of the utilized q = 1 PN-NR 
hybrid waveform. For higher mass-ratios, the PN-NR 
hybrids have a larger error due to the post-Newtonian 
waveform [21) . The error-bound on the hybrids increases 
with mass-ratio, however, is mitigated in our study here, 
because we use the q > 2 hybrids only for total mass of 
50 Mq , where less of the post-Newtonian waveform is in 
band. 

Because NR simulations are not available for arbitrary 
mass ratios, we will primarily concentrate our investiga- 
tion to the equal-mass and zero-spin NR hybrid wave- 
forms described above. The full IMR waveform can be 
generated at any point along the q — 1 line through a sim- 
ple rescaling of amplitude and phase with total mass M 
of the system. Despite such a simple rescaling, the q = 1 
line lies orthogonal to lines of constant chirp mass [22] , 
therefore tracing a steep gradient in terms of waveform 
overlap, and encompassing a large degree of waveform 
structure. 



C. PhenomB Waveforms 

Since our procedure for constructing an orthonormal 
basis begins with PhenomB waveforms, let us now in- 
vestigate how well these waveforms model the NR wave- 
forms to be interpolated. For this purpose, we adopt 
the notation hNp(M) and hpB(M) to represent NR and 
PhenomB waveforms of total mass M, respectively. We 
quantify the faithfulness of the PhenomB family by com- 
puting the mismatch M[h^(M) , 1ipb(M)] as a function 
of mass. The result of this calculation for 10 M < M < 
100 M Q is shown as the dashed curve in the top panel 



of Figure [T] The mismatch starts off rather high with 
Ai w 0.1 at 10 M and then slowly decreases as the mass 
is increased, until eventually flattening to M. s» 0.005 at 
high mass. 

The mismatch between NR and PhenomB waveforms 
can be reduced by optimizing over a mass-bias. This 
is accomplished by searching for the mass M' for which 
the mismatch A4[\inr(M), hps(M')] is a minimum. The 
result of this process is shown by the solid line in the top 
panel of Figure [l] Allowing for a mass bias significantly 
reduces the mismatch for M < 50 M . The mass M' that 
minimizes mismatch is generally smaller than the mass 
M of our NR "signal" waveform, M' < M over almost 
all of the mass range considered. Apparently, PhenomB 
waveforms are systematically underestimating the mass 
of the "true" NR waveforms, at least along the portion 
of parameter space considered here. The solid line in the 
bottom panel of Figure [T] plots the relative mass-bias, 
(M - M')/M. At 10 M this value is 0.3%, and it rises 
to just above 1% for 30 Mq. 1 

It is useful to consider how this mass bias compares 
to the potential parameter estimation accuracy in an 
early detection. For a signal with a matched-filter signal- 
to-noise ratio (SNR) of 8 — characteristic of early de- 
tection scenarios — template/waveform mismatches will 
influence parameter estimation when the mismatch is 
M > 8~ 2 / 2 ~ 0- 01 El- Placing a horizontal cut on 
the top panel of Figure [l] at M. = 8~ 2 /2, we see that for 
M > 40 M Q PhenomB waveform errors have no obser- 
vational consequence; for 15 M < M < 40 M a Phe- 
nomB waveform with the wrong mass will be the best 
match for the signal. For M < 15 M the missmatch 
between equal-mass PhenomB waveforms and NR (when 
optimizing over mass) grows to ~ 0.02. Optimization 
over mass-ratio will reduce this mismatch, but we have 
not investigated to what degree. 



III. INTERPOLATED WAVEFORM FAMILY 

A. PhenomB Template Bank 

We aim to construct an orthonormal basis via the SVD 
of a bank of PhenomB template waveforms, and then in- 
terpolate the coefficients of NR waveforms projected onto 
this basis to generate a waveform family with improved 
NR faithfulness. The first step is to construct a template 
bank of PhenomB waveforms, with attention restricted 
to equal-mass, zero-spin binaries. An advantage of focus- 
ing on the q = 1 line is that template bank construction 
can be simplified by systematically arranging templates 
in ascending order by total mass. 



In our calculation we have fixed q = 1 and x = 0. A more 
comprehensive minimization over mass, mass ratio, and effective 
spin might change this result. 
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With this arrangement we define a template bank to 
consist of N PhenomB waveforms, labelled g; = hpB(Mj) 
(i = 1,2, ... ,N), with M i+ i > M{ and with adjacent 
templates satisfying the relation: 



\0'-0( gi , gi+ i)\ <e, 



(5) 



where O' is the desired overlap between templates and 
e is some accepted tolerance in this value. The tem- 
plate bank is initiated by choosing a lower mass bound 
Mi = Af m in and assigning g! = h PB (Mi). Succes- 
sive templates are found by sequentially moving toward 
higher mass in order to find waveforms satisfying ([s]) un- 
til some maximum mass M max is reached. Throughout 
each trial, overlap between waveforms is maximized con- 
tinuously over phase shifts and discretely over time shifts. 
For template bank construction we choose to refine the 
optimization over time by considering shifts in integer 
multiples of At/100. 

We henceforth refer to our fiducial template bank 
which employs the parameters M m i n = 15 M©, M max = 
100 M , O' = 0.97, and e = 10~ 12 . The lower mass 
bound was chosen in order to obtain a reasonably sized 
template bank containing TV = 127 waveforms; pushing 
downward to 10 M© results in more than doubling the 
number of templates. Template waveforms each have a 
duration of 8 s and are uniformly sampled at At = s 
(a sample frequency of 32768 Hz). 508 MiB of memory 
is required to store this template bank using double- 
precision waveforms. 



B. Representation of Waveforms in a Reduced 
SVD Basis 

The next step is to transform the template waveforms 
into an orthonormal basis. Following the presentation in 
[S], this is achieved by arranging the templates into the 
rows of a matrix G and factoring through SVD to obtain 



G = VSU 1 , 



(6) 



where U and V are orthogonal matrices and S is a di- 
agonal matrix whose non-zero elements along the main 
diagonal are referred to as singular values. The SVD for 
G is uniquely defined as long as the singular values are 
arranged in descending order along the main diagonal of 
S. 

The end result of Q is to convert the N complex- 
valued templates into 27V real-valued orthonormal basis 
waveforms. The fc th basis waveform, uj., is stored in the 
k th row of U, and associated with this mode is the singu- 
lar value, <7/c, taken from the fc th element along the main 
diagonal of £. One of the appeals of SVD is that the 
singular values rank the basis waveforms with respect to 
their ability to represent the original templates. This can 
be exploited in order to construct a reduced basis that 
spans the space of template waveforms to some tolerated 
mismatch. 



For instance, suppose we choose to reduce the basis 
by considering only the first N' < 2N basis modes while 
discarding the rest. Template waveforms can be repre- 
sented in this reduced basis by expanding them as the 
sum 



JV' 

g' = ^2^ku k , 

k=l 



(7) 



where fik are the complex-valued projection coefficients, 

/i fe = (g,Ufc). (8) 

The prime in Q is used to stress that the reduced basis is 
generally unable to fully represent the original template. 2 
It was shown in [S] that the mismatch expected from 
reducing the basis in this way is 



(M) = (M(g',g)) 



1 

In 



2N 

E < 

k=N'+l 



(9) 



Given S, ([§ can be inverted to determine the number 
of basis waveforms, N', required to represent the original 
templates for some expected mismatch (M) . 

([9| provides a useful estimate to the mismatch in rep- 
reseting templates from a reduced SVD basis. In order to 
investigate its accuracy, however, we should compute the 
mismatch explicitly for each template waveform. Using 
the orthonormality condition (iij,Ufc) = Sjk, it is easy 
to show from Q that the mismatch between the tem- 
plate and its projection can be expressed in terms of the 
projection coefficients: 



X(g',g) = l- 



\ 



N' 

E 

fc=i 



(10) 



This quantity is minimized continuously over phase and 
discretely over time shifts in integer multiples of At. 

Choosing (M) = 10~ 6 , @ predicts that N' = 123 of 
the 2V = 254 basis waveforms from our fiducial tem- 
plate bank are required to represent the templates to the 
desired accuracy. In Figure [2] we compare the expected 
mismatch of 10 -6 to the actual mismatches computed 



from (10 1 for each PhenomB waveform in the template 
bank. The open squares in this plot show that the ac- 
tual template mismatch has a significant amount of scat- 
ter about (A4), but averaged over a whole remains well 
bounded to the expected result. The PhenomB template 
waveforms can thus be represented to a high degree from 
a reasonably reduced SVD basis. 

We are of course more interested in determining how 
well NR waveforms can be represented by the same re- 
duced basis of PhenomB waveforms. Since NR and Phe- 
nomB waveforms are not equivalent, ^ cannot be used 



2 In the case where TV' = 2N we are guaranteed from (j6j to com- 
pletely represent the template. 
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FIG. 2: Representation mismatch for the reduced SVD ba- 
sis with expected tolerance (M) = 10 -6 (traced by dot- 
ted line). Open squares (open circles) show the represen- 
tation mismatch for PhenomB (NR) waveforms evaluated at 
the PhenomB template locations, while the thin solid line 
(thin dashed line) traces representation mismatch of Phe- 
nomB (NR) waveforms evaluated between templates. The 
NR waveforms cannot be represented as well as their Phe- 
nomB counterparts, although their total match is improved 
over using PhenomB waveforms alone. This is evidenced by 
the thick solid line tracing the mass-optimized mismatch be- 
tween NR and PhenomB waveforms (i.e. the solid line in the 
top panel of Figure fill. 
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FIG. 3: Convergence of representation mismatch with tight- 
ening tolerance (M). As a function of (M), we plot the 
averages of the four data-sets shown in Figure [2] The rep- 
resentation mismatch of the PhenomB waveforms for masses 
in the template bank used to construct the SVD basis decays 
roughly with (M). The representation mismatch of PhenomB 
waveforms for masses between the masses in the template bank 
reaches a plateau of ~ 10 -5 (the precise value depends on 0', 
cf. |5|). The representation mismatch of NR waveforms is 
yet larger with a plateau of ~ 10~ 3 that flattens for larger 
(M) (this flattening depends only mildly on O'). 



to estimate the mismatch obtained when projecting NR 
waveforms onto the reduced basis. We must therefore 
compute their representation mismatch explicitly. A gen- 
eral waveform, h, can be represented by the reduced basis 
in analogy to |7| by expressing it as the sum: 



N' 



(11) 



k=l 



where fit- = (h, Ufc) . As before, the represented waveform 
h' will in general be neither normalized nor equivalent to 
the original waveform h. The mismatch between them is 



M(h',h) = 1 



\ 



N' 



X] wit, 



(12) 



fc=i 



where we remind the reader that we always minimize over 
continuous phase shifts and discrete time shifts of the two 
waveforms. 

In Figure [2] we use open circles to plot the representa- 
tion mismatch of NR waveforms evaluated at the same 
set of masses Mj from which the PhenomB template bank 
was constructed. We see that NR waveforms can be rep- 
resented in the reduced basis with a mismatch less than 
10~ 3 over most of the template bank boundary. This 
is about a factor of five improvement in what can be 
achieved by using PhenomB waveforms optimized over 
mass. Since NR waveforms were not originally included 
in the template bank, and because a mass-bias exists 



between the PhenomB waveforms which were included, 
we can expect that the template locations have no spe- 
cial meaning to NR waveforms. This is evident from 
the thin dashed line which traces the NR representation 
mismatch for masses evaluated between the discrete tem- 
plates. This line varies smoothly across the considered 
mass range and exhibits no special features at the tem- 
plate locations. This is in contrast to the thin solid line 
which traces PhenomB representation mismatch evalu- 
ated between templates. In this case, mismatch rises as 
we move away from one template and subsequently falls 
back down as the next template is approached. 

The representation tolerance (M) of the SVD is a 
free parameter, which so far, we have constrained to 
be (Ai) = 10~ 6 . When this tolerance is varied, we ob- 
serve the following trends: (i) PhenomB representation 
mismatch generally follows (M); (ii) NR representation 
mismatch follows (M) at first and then saturates to a 
minimum as the representation tolerance is continually 
reduced. These trends are observed in Figure [3] where 
we plot NR and PhenomB representation mismatch av- 
eraged over the mass boundary of the template bank 
evaluated both at and between templates. The satura- 
tion in NR representation mismatch occurs when the re- 
duced basis captures all of the NR waveform structure 
contained within the PhenomB basis. Reducing the ba- 
sis further hits a point of diminishing returns as the in- 
creased computational cost associated with a larger basis 
outweighs the benefit of marginally improving NR match. 
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C. Interpolation of NR Projection Coefficients 

We now wish to examine the possibility of using the 
reduced SVD basis of PhenomB template waveforms to 
construct a new waveform family with improved NR rep- 
resentation. The new waveform family would be given by 
a numerical interpolation of the projection coefficients of 
NR waveforms expanded onto the reduced basis. Here 
we test this using the fiducial template bank and reduced 
basis described above. 

The approach is to sample NR projection coefficients, 
fj,k(xi), at some set of locations, Xi, and then perform 
an interpolation to obtain the continuous function fj/ k (x) 
that can be evaluated for arbitrary x. The accuracy 
of the interpolation scheme is maximized by finding the 
space for which Hk{ x i) are smooth functions of x. It is 
reasonable to suppose that the projection coefficients will 
vary on a similar scale over which the waveforms them- 
selves vary. Hence, a suitable space to sample along is 
the space of constant waveform overlap. We define this to 
be the space x — [—1,1] for which the physical template 
masses are mapped according to: 



Mi -^ Xi = -l 



t i-l 

'N -1 



(13) 



Moving a distance Ax = 2/(N — 1) in this space is thus 
equivalent to moving a distance equal to the overlap be- 
tween adjacent templates. 

In this space, we find the real and imaginary compo- 
nents, lZ/ik(x) and X/jj c (x), of the complex projection 
coefficients to be oscillatory functions that can roughly 
be described by a single frequency. This behaviour is 
plotted for the basis modes k = 1, 50, and 123 in Fig- 
ure [4] Another trend observed in this plot is that the 
projection coefficients become increasingly complex (i.e. 
show less structure) for higher-order modes. This is a 
direct result of the increasing complexity of higher-order 
basis waveforms themselves. We find that the low-order 
waveforms are smoothest while the high-order modes fea- 
ture many of the irregularities associated with the mul- 
tiple frequency components and merger features of the 
templates. Though they are more complex, higher-order 
modes have smaller singular values and are therefore less 
important in representing waveforms. This is evident 
from the steady decline in amplitude of the projection 
coefficients when moving down the different panels of Fig- 
ure m 

We shall use Chebyshev polynomials to interpolate the 
projection coefficients. These are a set of orthogonal 
functions where the j th Chebyshev polynomial is defined 
as 



Tj(x) = cos(j arccosz), 



e[-l,l] 



(14) 



The orthogonality of Chebyshev polynomials can be ex- 
ploited to perform an n th order Chebyshev interpolation 
by sampling Hk{x) at the n+1 so-called collocation points 
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FIG. 4: Plotted are the real part of the complex coefficients 
obtained from projecting NR waveforms onto the PhenomB 
SVD basis waveforms ui (top panel), U50 (third panel), and 
Ui23_(fifth panel). The x axis has been constructed according 



to ( 13 1 and is ideal for performing a Chebyshev interpolation. 



Open circles show the projection coefficients sampled at the 
collocation points in ( |15[ ) for ann — 175 Chebyshev interpola- 
tion; solid lines trace the resultant interpolation for each basis 
mode. Interpolation is performed separately on the real and 
imaginary parts of y,k{x) and combined afterward to obtain 
the complex function /Xfe(x). Below each set of coefficients is 
the absolute value of the interpolation residual \r^\ defined in 
pi. 



given by the Gauss-Lobatto Chebyshev nodes [24] 



-cos 



(15) 



for i = 0, 1, . . . , n. 

In general, the interpolation will not be exact and some 
residual, rfc(M), will be introduced: 



r fe (M)=/4(M)- Mfe (M). 



(16) 



Here /^(Af) is the actual coefficient of 1inr(M) projected 
onto the basis waveform u,t, 



ft k (M) = (h NR (M),u fe ), 



(17) 



and n' k (M) is the coefficient obtained after interpolation. 
The new waveform family is expressed numerically as a 
function of mass through the relation 



N' 



N 



htatp(M) = ^'k( M ) u k = i^(M) + r k {M)\ u fc) 

(18) 



k=l 



k=l 
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where the subscript "intp" reminds the reader that this 
is computed from an interpolation over An inter- 
polated waveform of total mass M can be compared to 
the original NR waveform (which we consider to be the 
"true" signal), where the latter is expressed as 



inr 



(M) 



2N 

E 

k=l 



H k (M)u k + h ± (M), 



(19) 



with h±(M) denoting the component of Iinr(M) that is 
orthogonal to the SVD basis (i.e. orthogonal to all Phc- 
nomB waveforms in the template bank). Iinr(M) differs 
from hi n tp(-M) by an amount 



Sh = h intp (M) - h NR (M) 



N 

E 



r fc (M)u fe 



2A r 

■ E 

fc=jV' + l 



Hk(M)u k -h±{M) (20) 



To compute the impact of the various approxima- 
tions influencing (201, we calculate the overlap between 



the interpolated waveform, and the exact waveform, 
O [h; ntp (M), hNB,(M)]. To begin this calculation, it is 
useful to consider the square of the overlap, 



0(h + <5h,h) 



(h,h + (5h)(h + (5h,h) 
(h,h)(h + <5h,h + <5h)' 



(21) 



where we have dropped the explicit mass-dependence and 
subscripts for convenience. Using (h, h) = 1 and Taylor- 



expanding the right-hand-side of (21 ) to second order in 
Sh, we find 

0(h + Sh, h) 2 = 1 - (Sh, Sh) + (Sh, h) (h, Sh) . (22) 
To second order in Sh, the mismatch is therefore 

M(h + Sh, h) = hsh, Sh) - I (Sh, h) (h, Sh) . (23) 



We note that the right-hand-side of ( 23 ) can be written 
as |(<5hj_, <5h_i_), where Sh± is the part of Sh orthogonal 
to h, 

Shx =5h-^h,h)h. (24) 
However, for simplicity, we proceed by dropping the last 



term in ( 23 1 



M(h + 5h,h) < -{6h,5h). 



(25) 



Using (20), this gives 



X[h intp (M),h NR (M)] < 



1 N ' 

oEM M )l 

k=l 



2JV 

E 

fc=jV' + l 



(26) 



We thus see three contributions to the total mismatch: 
(i) the interpolation error, X^feLi l r fc(-^0| 2 ; (ii) the trun- 
cation error from the discarded waveforms of the reduced 




120 



FIG. 5: The solid line traces the interpolation error R\ in ( 27 1 
maximized over mass for each mode k. This is used in (281 
to define an upper bound to the interpolation error arising in 
(261. The cumulative sum in the latter expression is traced 
by the dashed curve and shows that the interpolation error is 
dominated by the lowest-order basis waveforms. 



basis, J2k=N'+i \^k{M)\ 2 ; (iii) the failure of the SVD ba- 
sis to represent the NR waveform, |hjj. The sum of the 
last two terms, which together make up the representa- 
tion error, is traced by the dashed line in Figure [2j The 
goal for our new waveform family is to have an interpo- 
lation error that is negligible compared to the represen- 
tation error. 

To remove the mass-dependence of interpolation error 
in ( |2G| , we introduce the maximum interpolation error of 
each mode, 



R k = max | r k {M) | . 



This allows the bound 



N' N' 



(27) 



(28) 



k=l 



k=l 



to place an upper limit on the error introduced by inter- 
polation. Figure[5]plots R\ as a function of mode-number 

k as well as the cumulative sum 

J2k=i R k/ 2 - The data 
pertains to an interpolation performed using n = 175 

Chebyshev polynomials on the reduced SVD basis con- 
taining the frist N' — 123 of 2N — 254 waveforms. In 
this case, we find the interpolation error to be largely 
dominated by the lowest-order modes and also partially 
by the highest-order modes. Interpolated coefficients for 
various modes are plotted in Figure [4] and help to explain 
the features seen in Figure [5j In the first place, interpola- 
tion becomes increasingly more difficult for higher-order 
modes due to their increasing complexity. This prob- 
lem is mitigated by the fact that high-order modes are 
less important for representing waveforms, as evidenced 
by the diminishing amplitude of projection coefficients. 
Although low-order modes are much smoother and thus 
easier to interpolate, their amplitudes are considerably 
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FIG. 6: Open circles show the total mismatch of our new 
waveform family obtained by interpolating the coefficients of 
NR waveforms projected onto a reduced SVD basis of Phe- 
nomB waveforms. To highlight the error introduced by in- 
terpolation, the solid curve traces only the mismatch of NR 
waveforms represented by the reduced basis (i.e. the dashed 
curve in Figure [2b. The total interpolation mismatch is lower 
than the dotted line tracing the mass-optimized mismatch be- 
tween NR and PhenomB waveforms (i.e. the solid line in the 
top panel of Figure [I]) and demonstrates the ability of SVD 
to boost the accuracy of the PhenomB family. 



larger meaning that interpolation errors are amplified 
with respect to high-order modes. 



( 26 ) summarizes the three components adding to the fi- 



nal mismatch of our interpolated waveform family. Their 
total contribution can be computed directly from the in- 
terpolated coefficients in a manner similar to (12): 



A*[h intp (M),h NR (M)] = 1- 



\ 



N' 



E^WW' ( 29 ) 



k = l 



In the case of perfect interpolation for which /jf k (M) — 
jUfc(M), ([18| and @ reduce to (fill and Q respec- 



tively, and the total mismatch is simply the representa- 
tion error of the reduced basis. 

In Figure [6] open circles show the total mismatch ( 29 ) 
between our interpolated waveform family and the true 
NR waveforms for various masses. Also plotted is the 
NR representation error without interpolation and the 
mismatch between NR and PhenomB waveforms min- 
imized over mass. We see that interpolation introduces 
only small additional mismatch to the interpolated wave- 
form family, and remains well below the optimized NR- 
PhenomB mismatch. This demonstrates the efficacy of 
using SVD coupled to NR waveforms to generate a faith- 
ful waveform family with improved accuracy over the ef- 
fectual PhenomB family that was originally used to cre- 
ate templates. This represents a general scheme for im- 
proving phenomenological models and presents an inter- 
esting new opportunity to enhance the matched-filtering 
process employed by LIGO. 



IV. HIGHER DIMENSIONS 

So far, we have focused on the total mass axis of pa- 
rameter space. As already discussed, this served as a 
convenient model-problem, because the q = 1 NR wave- 
form can be rescaled to any total mass, so that we are 
able to compare against the "correct" answer. The nat- 
ural extension of this work is to expand into higher di- 
mensions where NR waveforms are available only at cer- 
tain, discrete mass-ratios q. In this section we consider 
expanding our approach of interpolating NR projection 
coefficients from a two-dimensional template bank con- 
taining unequal-mass waveforms. 

We compute a template bank of PhenomB waveforms 
covering mass-ratios q from 1 to 6 and total masses 
50 M Q < M < 70 M Q . This mass range is chosen to fa- 
cilitate comparison with previous work done by |12j . For 
the two-dimensional case the construction of a template 
bank is no longer as straightforward as before due to the 
additional degree of freedom associated with varying q. 
One method that has been advanced for this purpose is 
to place templates hexagonally on the waveform mani- 
fold [25]. Using this procedure we find N = 16 templates 
are required to satisfy a minimal match of 0.97. 

Following the waveform preparation of |12j , templates 
are placed in the rows of a matrix G with real and imag- 
inary components filled in alternating fashion where the 
whitened waveforms are arranged in such a way that 
their peak amplitudes are aligned. The waveforms are 
sampled for a total duration of 2 s with uniform spac- 
ing At = 2~ 15 s so that 16 MiB of memory is required to 
store the contents of G if double precision is desired. Ap- 
plication of (|6| transforms the 16 complex-valued wave- 
forms into 32 real- valued orthonormal basis waveforms. 

The aim is to sample the coefficients of NR wave- 
forms projected onto the SVD basis of PhenomB wave- 
forms using mass-ratios for which NR data exists, and 
then interpolate amongst these to construct a numeri- 
cal waveform family that can be evaluated for arbitrary 
parameters. This provides a method for evaluating full 
IMR waveforms for mass-ratios that have presently not 
been simulated. To summarize, we take some NR wave- 
form, hiMfjXM, q), or total mass M and mass-ratio q, and 
project it onto the basis waveform Ufe in order to obtain 



fi>k(M,q) = (h NR (M,g),u fc ). 



(30) 



Next we apply some two-dimensional interpolation 
scheme on (30) to construct continuous functions 



fj,' k (M, q) that can be evaluated for arbitrary values of 
M and q bounded by the regions of the template bank. 
The interpolated waveform family is given numerically 
by the form: 



N' 



h in tp(M,g) = ^2(i k (M,q)u k . 



(31) 



k=l 



As before, the interpolation process works best if we 
can develop a scheme for which the projection coefficients 
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FIG. 7: (left panels) Real components of the smoothed PhenomB projection coefficients fbz{M, q) (top) and fi^^M, q) (bottom) 
for mass-ratios 1 < q < 6 and total masses 50 Mq < M < 70 Mq. (middle panels) Real components of the smoothed NR 
projection coefficients jls{M,q) (top) and /Etie(M, q) (bottom) for mass-ratios q = {1, 2, 3, 4, 6} and total mass 50 Mq < 
M < 70 Mq. (right panels) Real components of the smoothed PhenomB projection coefficients jl^{M,q) (top) and flie(M,q) 
(bottom) coarsened to the same set of mass-ratios as the middle panels. For both the middle and rightmost panels an artificial 
row of white space has been plotted for q — 5 in order to ease comparison with the leftmost panels. 



are smoothly varying functions of M and q. Following 
the procedure described in [12] , the complex phase of the 
first mode is subtracted from all modes: 



(32) 



To motivate why (32) might be useful, let us con- 



sider modifying the PhenomB waveform family with a 
parameter-dependent complex phase $(M, q): 



h PB (M,q) ^ e^ M ^hp B (M,q). 



(33) 



When constructing a template bank, or when using a 
template bank, such a complex phase $(M, q) is irrele- 
vant, because the waveforms are always optimized over a 
phase-shift. However, $(M, q) will appear in the projec- 
tion coefficients, (30), 



Vk(M, q) — »• e i9 ^ M,q ^fj,k(M,q). 



(34) 



Therefore, if one had chosen a function $(M, q) with 
fine-scale structure, this structure would be inherited by 
the projection coefficients /ifc (M, q) . For traditional uses 
of waveform families the overall complex phase $(M, q) 
is irrelevant, and therefore, little attention may have 
been paid to how it varies with parameters (M, q) . The 
transformation ( 32 ) removes the ambiguity inherent in 
$(M, q) by choosing it such that arg/ii(M, q) = 0. This 
choice ties the complex phase to the physical variations 
of the /ii coefficient, and does therefore eliminate all un- 
physical phase- variations on finer scales. 

In the leftmost panels of Figure [7] we plot the real part 
of the smoothed coefficients jlk (M, q) for PhenomB wave- 
forms projected onto the basis modes k — 3 and k = 16. 
The middle panels show the same thing except using the 
NR waveforms evaluated at the set of mass-ratios q — 



{1, 2, 3, 4, 6} for which we have simulated waveforms. 
Obviously, the refinement along the q axis is much finer 
for the PhenomB waveforms since they can be evaluated 
for arbitrary mass-ratio, whereas we arc limited to sam- 
pling at only 5 discrete mass-ratios for NR waveforms. 
For comparison purposes, the rightmost panels of Figure 
[7] show the PhenomB projection coefficients coarsened to 
the same set of mass-ratios for which the NR waveforms 
are restricted to. 



We find the same general behaviour as before that 
low-order modes display the smoothest structure, while 
high-order modes exhibit increasing complexity. A plau- 
sible interpolation scheme would be to sample jlk for 
NR waveforms of varying mass for constant mass ratio 
(i.e. as we have done previously) and then stitch these 
together across the q axis. Since the projection coeffi- 
cients in Figure [7] show sinusoidal structure they must be 
sampled with at least the Nyquist frequency along both 
axes. However, looking at the middle and rightmost pan- 
els it appears as though this is not yet possible given the 
present set of limited NR waveforms. At best the 5 avail- 
able mass-ratios are just able to sample at the Nyquist 
frequency along the q axis for high-order modes. In order 
to achieve a reasonable interpolation from these projec- 
tion coefficients the current NR data thus needs to be 
appended with more mass-ratios. Based on the left pan- 
els of Figure [7] a suitable choice would be to double the 
current number of mass-ratios to include q = {1.5, 2.5, 
3.5, 4.5, 5, 5.5}. Hence, though it is not yet practical to 
generate an interpolated waveform family using the SVD 
boosting scheme applied to NR waveforms, the possibil- 
ity remains open as more NR waveforms are generated. 
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V. DISCUSSION 

We have shown that SVD can be used to improve the 
representation of NR waveforms from a PhenomB tem- 
plate bank. A reasonably reduced SVD basis was able 
to reduce mismatch by a factor of five compared to Phe- 
nomB waveforms optimized over mass. There was also no 
mass-bias associated with the SVD basis and therefore no 
optimization over physical parameters required. This oc- 
curs because SVD unifies a range of waveform structure 
over an extended region of parameter space so that any 
biases become blended into its basis. SVD therefore rep- 
resents a generalized scheme through which phenomeno- 
logical waveform families can be de-biased and enhanced 
for use as matched-filter templates. 

We were able to calibrate an SVD basis of Phe- 
nomB templates against NR waveforms in order to con- 
struct a new waveform family with improved accuracy. 
This was accomplished by interpolating the coefficients 
of NR waveforms projected onto the PhenomB basis. 
Only marginal error was introduced by the interpola- 
tion scheme and the new waveform family provided a 
more faithful representation of the "true" NR signal com- 
pared to the original PhenomB model. This was shown 



explicitly for the case of equal-mass, zero-spin binaries. 
We proceeded to investigate the possibility of extend- 
ing this approach to PhenomB template banks contain- 
ing unequal-mass waveforms. At present, however, this 
method is not yet feasible since the current number of 
mass-ratios covered by NR simulations are unable to sam- 
ple the projection coefficients with the Nyquist frequency. 
This method will improve as more NR waveforms are sim- 
ulated and should be sufficient if the current sampling 
rate of mass-ratios were to double. 
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